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Chaos synchronization in a ring of diffusively coupled nonlinear oscillators driven by an exter- 
nal identical oscillator is studied. Based on numerical simulations we show that by introducing 
additional couplings at (mNc + l)-th oscillators in the ring, where m is an integer and Nc is the 
maximum number of synchronized oscillators in the ring with a single coupling, the maximum num- 
ber of oscillators that can be synchronized can be increased considerably beyond the limit restricted 
by size instability. We also demonstrate that there exists an exponential relation between the num- 
ber of oscillators that can support stable synchronization in the ring with the external drive and 
the critical coupling strength Ec with a scaling exponent 7. The critical coupling strength is calcu- 
lated by numerically estimating the synchronization error and is also confirmed from the conditional 
Lyapunov exponents (CLEs) of the coupled systems. We find that the same scaling relation exists 
for m couplings between the drive and the ring. Further, we have examined the robustness of the 
synchronous states against Gaussian white noise and found that the synchronization error exhibits a 
power-law decay as a function of the noise intensity indicating the existence of both noise-enhanced 
and noise-induced synchronizations depending on the value of the coupling strength e. In addi- 
tion, we have found that Sc shows an exponential decay as a function of the number of additional 
couplings. These results are demonstrated using the paradigmatic models of Rossler and Lorenz 
oscillators. 

PACS numbers: 05.45.Xt 



I. INTRODUCTION 

Chaos synchronization has been receiving a great deal 
of interest for more than three decades [iH5| . In partic- 
ular, chaos synchronization in arrays of coupled nonlin- 
ear dynamical systems has been extensively investigated 
over the years in view of its diverse applications in spa- 
tially extended systems, neural process, networks, etc. 
[6l4ll|. Linear arrays with periodic boundary condition 
(ring geometry) have been used widely in modeli ng phys - 



iological, biochemical and biological phenomena 111 -[T3|. 
For example, morphogenesis in biological context |14| and 
transitions between different animal gaits have been ex- 
plained by considering a model composed of a ring of cou- 
pled oscillators [l^ . An important application of the ring 
geometry is that the resulting spatiotemporal patterns 
in the ensemble of coupled oscillators can be analyzed 
through symmetry arguments [T^ . ITsj . Recently several 
interesting dynamical properties/collective behaviors in- 
cluding amplitude death and chimera states have been 
identified in such a ring type configuration [TT| - [l3l HgI - 

Some of the recent studies have considered synchro- 
nization dynamics in both ring and linear arrays cou- 
pled together in order to understand the dynamics of 
basic units of networks [111, [l8l - l20j . Recently, interest- 
ing scaling behavior of correlation properties of interact- 
ing dynamical s;Y^tems in such a configuration has been 
demonstrated [20[. However, most of the studies have 
considered unidirectional coupling in both the ring and 
linear arrays. Because of the diverse nature of interaction 



in real world phenomena, we have considered diffusively 
(nearest-neighbor) coupled chaotic systems with ring ge- 
ometry driven by an external identical drive in view of 
its widespread applications in engineering, robotics, net- 
works, and physiological and biological systems (25l. [26|. 
For instance, cultured networks of heart cells are exam- 
ples of bi olog ical structures with strong nearest-neighbor 
couphng [11 113. 



The phenomenon of size instability, where a critical 
size of the number of oscillators upto which a stable syn- 
chronous chaotic state exists, of a uniform synchronous 
state in arrays of coupled oscillators with both periodic 
and free-end boundary conditions have been widely stud- 
ied @, 0, [13: 0; [28l - [33j . Increasing the number of oscil- 
lators beyond this limit leads to desynchronization and 
the occurrence of spatially incoherent behavior (eg., high- 
dimensional or spatiotemporal chaos). The stability of 
synchronous chaos in coupled dynamical systems plays a 
crucial role in the study of p attern formation, spatiotem- 
poral chaos, etc. 0, Isol - fsa [13, [s^. In this connection, 
using the paradigmatic models of Rossler and Lorenz os- 
cillators wc shall demonstrate in this paper that the max- 
imal number of oscillators in the ring geometry that can 
support stable synchronous chaos can be increased by 
integer multiples of the original number of oscillators in 
the ring with additional couplings from the same drive 
oscillator. Furthermore, it is also found that the critical 
coupling strength and the number of oscillators which 
can be synchronized in the ring exhibit an exponential 
relation with a scaling exponent and indeed the relation 
remains unaltered even on increasing the number of cou- 
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plings between the drive and the ring. In addition, the 
synchronization error displays a power law decay as a 
function of the noise intensity for a fixed value of the 
coupling strength indicating the existence of noise en- 
hanced/induced synchronization. It is to be noted that 
small world networks can be generated by introducing 
additional couplings between randomly selected nodes to 
create shortest paths (links) between distant nodes [lo| . 
Further, recent studies on synchroniz ability of networks 
have been employing pinning control in which hubs in 
the networks are connected to the same drive node [l^l 
as in our present study. 

In particular, we consider here the Rossler and Lorenz 
oscillators in a ring geometry with diffusive coupling be- 
tween them and driven by an external identical oscil- 
lator, whose strength is proportional to a parameter e 
(see Eq. (2) below). Based on numerical simulations, 
we find that the critical coupling strength, say Ec, below 
which no synchronization exists (e < Ec), of the exter- 
nal drive increases exponentially with a scaling exponent, 
7 € (0.3,0.5), as a function of the number of oscillators 
in the array that supports a stable synchronous state. 
Further we observe that the number of oscillators which 
supports such a stable synchronous state can be increased 
in integer multiples by introducing additional couplings 
at the {mNc + l)-th oscillator of the ring with same value 
of the coupling strength, where m is the number of cou- 
plings and Nc is the maximum number of oscillators in 
the ring that can sustain stable synchronization with a 
single coupling. Interestingly, this exponential relation is 
maintained while increasing the number of couplings, m, 
between the array and the external drive. In addition, 
we have found that £c shows an exponential decay as a 
function of the number of additional couplings between 
the drive and the response array for a fixed number of os- 
cillators in the array. Further, we find that these results 
are robust against Gaussian white noise of small intensity 
and the synchronization error exhibits a power law decay 
as a function of the noise intensity. These results also 
indicate the existence of both the phenomena of noise- 
enhanced and noise-induced synchronizations depending 
on the value of the coupling strength e beyond certain 
threshold values of the noise intensity. It is to be empha- 
sized that all the numerical simulations throughout the 
manuscript have been repeated with several initial con- 
ditions and the results are the averages of a large number 
of realizations. 

The structure of the paper is as follows. In Sec. |IT1 
we study synchronization in rings of diffusively coupled 
Rossler and Lorenz systems driven by external identical 
oscillators with a drive-response configuration. We show 
that in both the cases the critical coupling strength in- 
creases exponentially with the number of oscillators in 
the response array with a scaling exponent. In addition, 
the systems exhibit a power law decay of the synchro- 
nization error as a function of the noise intensity demon- 
strating noise-enhanced and noise-induced synchroniza- 
tions depending on the value of the coupling strength. 



In Sec. mil we demonstrate that the size of the ring can 
be increased beyond the size instability limit by integer 
multiples of the maximum number of synchronized oscil- 
lators (Nc) in the ring with a single coupling for the same 
value of coupling strength. We also show that the same 
scaling relation with a characteristic exponent is valid for 
the case of arbitrary number of couplings also. Further, 
we examine the effect of noise on the robustness of the 
synchronous state for a fixed value of the noise intensity 
and as a function of the noise intensity for a fixed value of 
the coupling strength in all the cases. Finally, in Sec. lIVi 
we present a summary and conclusions. 



II. CHAOS SYNCHRONIZATION IN 
DIFFUSIVELY COUPLED OSCILLATORS 
DRIVEN BY AN EXTERNAL IDENTICAL 
OSCILLATOR 

For the present study, we consider a coupling scheme 
with drive-response configuration as shown in Fig. [U in 
which the diffusively coupled circular array is driven by 
an external identical oscillator. Here '0' denotes the ex- 
ternal drive and '1', '2', . . . W denote the constituents of 
the response array with nearest neighbor (diffusive) cou- 
pling. For simplicity, we assume that all the oscillators 
in the drive-response configuration are identical. In this 
section, we study chaos synchronization in the response 
array when driven by the external drive. For this pur- 
pose, we have considered Rossler and Lorenz oscillators 
coupled according to Fig. [TJ The value of the diffusive 




FIG. 1: (color online) Schematic diagram of the ring of diffu- 
sively coupled oscillators driven by an external identical os- 
cillator with drive-response configuration. 

coupling constant is chosen such that all the oscillators 
in the array are in a stable synchronous state. By vary- 
ing the number of oscillators in the response array from 
N = 2 onwards, we calculate the critical value of the cou- 
pling constant Ec at and above which the response array 
evolves in synchrony with the external drive. 
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A. Coupled Rossler system 

First we analyse the coupled Rossler systems given by 
the dynamical equations 



io= - (yo + zo), 
zo=b + zn{xn - c), 



(1) 



+ d{xj+i + x-j^x — 2xj) + 5ije{xo — x-j), (2a) 
Vj = + m (2b) 
z,^b + zj{xj-c), J - 1,2,--- ,iV, (2c) 



where a = 0.15, b = 
the coupling strength. 



0.2, c = 10, d = 1 and e is 
Here the variable xq, yo and 
zq correspond to the drive system and Xj, yj and Zj 
(j = 1,2,...,A^) represent the diffusively coupled re- 
sponse array. The first oscillator of the response array 



([2]) is driven by the drive (P). 
function given by 



Si , is the Kronecker delta 



1 for i = j 
otherwise. 



The isolated system ([T]) exhibits chaotic behavior for the 
above choice of parameters with the Lyapunov exponents 
Ai « 0.1304 > 0, As = and A3 w -14.1405. 

Next we shall study the dynamics of the drive-response 
configuration ([T]) and . The value of the diffusive cou- 
pling constant d is chosen such that it supports the maxi- 
mum number of oscillators in the circular array to evolve 
in synchronous fashion. In the absence of any external 
coupling (e = 0), the array of diffusively coupled Rossler 
systems ([2|) exhibits chaos synchronization in all the N 
oscillators with N < iVmax with Q 



(\/Amax/4d 



sm 



(3) 



where Amax is the maximal Lyapunov exponent. One 
can easily check that iVniax = 17 for the above choice of 
parameters. 

As soon as the external drive is switched on (e ^ 0), 
the stable synchronous state of the diffusively coupled os- 
cillators gets destroyed for small values of e. The number 
of oscillators in the response array ([2]) which retains syn- 
chronization with the drive ([1} depends on the coupling 
strength e. To estimate the quality of synchronization, 
we define a measure, namely the synchronization error, 
as the Euclidean norm 



1 

N 



N 



1 
2 

(4) 



In order to get a perfect synchronization in the drive- 
response configuration, ([1]) and we require 77 as 
t — !■ 00. We remark here that all the simulations in the 
manuscript are performed for an average of over 100 ran- 
dom initial conditions. By examining 77, we extract the 
critical coupling strength £c corresponding to the number 
of oscillators, N , in the response array, whose dynam- 
ics arc entrained with that of the drive in phase space. 
For example, for — 2, we find the critical coupling 
strength as = 0.47. The £c values are given in Table U 
for different values of N . However, for N > Q, the system 
gets destabilized or becomes completely unstable, for any 
value of e 7^ 0. 

TABLE I: Critical coupling strength Sc 



Rossler system 


Lorenz system 


iV Ec 


iV 




2 0.447 


2 


12.528 


3 0.721 


3 


20.288 


4 1.134 


4 


31.334 


5 1.882 


5 


51.486 



By careful examination of the numerical values, we re- 
alize that there is an exponential relation between the 
number of oscillators (A^) in the response array that sup- 
ports synchronization with the external drive and the 
critical coupling strength, ec- By fitting the numerical 
values, one can easily establish the relation connecting ec 
and A^ as 



£oexp(7A^), 



(5) 



with the proportionality constant £0 = 0.1665 and the 
scaling exponent 7 = 0.4842. The critical coupling 
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FIG. 2: (color online) Variation of the critical coupling 
strength, Ec, as a function of the number of oscillators, A'^, 
for (a) Rossler system (eo ~ 0.1665 and 7 = 0.4842) and (b) 
Lorenz system {eq = 4.8158 and 7 = 0.4731). Filled circles 
correspond to the numerical data and solid lines are the plot 
of Ec using relation ((5)1. 

strength as a function of the number of oscillators is 
shown in Fig. [2{a), where the filled circles correspond 
to numerical data and the solid line is the fitted data us- 
ing the relation ([S]) . The above result is also confirmed by 
calculating the conditional Lyapunov exponents (CLEs) 
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maximum of iV = 9 oscillators in the stable synchronous 
state as per Eq. 

When the external coupling is switched on, the syn- 
chronization in the array (O gets lost. However, by 
choosing the number of oscillators, iV, and the coupling 
strength, e, appropriately one can synchronize the array 
with the external drive. The values of Ec for different 
numbers of oscillators in the response array, N, arc also 



FIG. 3: (color online) Conditional Lyapunov exponents of the 
coupled Rossler system ((2)) for different values of the number 
of oscillators (A'^) in the array as a function of the coupling 
strength. 



of the response array by increasing the number of the 
coupled oscillators N. Here by CLE we mean the largest 
Lyapunov exponent of the response array, which of course 
should be less than zero in order to have a synchronous 
state. Fig. [3] shows the variation of CLEs as a function 
of £ for the various numbers of oscillators in the response 
array. From Fig. |31 one can identify that the value of Ec 
for each A'', where a transition of the CLE from positive 
to negative value occurs, agrees with the value of Ec given 
in Table L 



B. Coupled Lorenz systems 



Next we consider the case of the Lorenz 
scribed by the following set of equations [36 
drive, 

Xq =(T{yo - xq), 

yo = rxo -yo- XqZo 

zo = - l3zo + Xoyo, 



stem dc- 
as the 



(6) 



where tT = 10, r = 23 and /3 = 1. A ring of diffusively 
coupled Lorenz systems for the response array can be 
represented by the following set of coupled equations. 



- + Xj^i — 2xj 



Pzj +Xjyj, j = 1,2,--- ,iV, 



(7a) 
(7b) 
(7c) 



where d is the diffusive coupling constant and e is the 
coupling strength of the external drive. One may note 
that the isolated system ^ exhibits chaotic behavior for 
the above choice of parameters with the Lyapunov expo- 
nents Ai 0.6075 > 0, A3 and A3 = -17.9194. 

We study the variation of the critical coupling strength 
as a function of the number of oscillators in the response 
array by numerically examining the synchronization error 
77 in a similar fashion as in the case of the coupled Rossler 
oscillators. For this purpose, we fix the diffusive coupling 
constant as d = 30. In this case, the response array 
([7]), in the absence of external drive (e = 0), supports a 




FIG. 4: (color online) Conditional Lyapunov exponents of the 
coupled Lorenz system ((7| for different values of the number of 
oscillators in the array as a function of the coupling strength. 

given in Table HI Fig. O^b) depicts the plot of TV versus 
Ec- It is easy to see that, as in the case of the coupled 
Rossler oscillators discussed above, Ec again increases ex- 
ponentially as a function of N according to the relation 
(P with £0 = 4.8158 and 7 = 0.4731. 

The above results have also been examined by calculat- 
ing the conditional Lyapunov exponents associated with 
d?]) and it is confirmed that the largest CLE (cf. Fig. |4]) 
transits from positive to negative value at the critical 
coupling strength of a given N as shown in Table HI 

So far. we have considered Rossler and Lorenz oscilla- 
tors with a ring geometry and driven by an external iden- 
tical oscillator in a drive-response configuration. Based 
on numerical simulations, we have found that there is an 
exponential relation connecting the critical value of the 
coupling strength Ec and the number of oscillators in the 
response array N that evolve in synchrony with the ex- 
ternal drive with a scaling exponent 7 0.48 for both 
the Rossler and Lorenz oscillators. Next, we will exam- 
ine the effect of Gaussian white noise on the stability of 
synchronization of the coupled oscillators. 



C. Effect of noise 

We have examined the robustness of the synchronous 
states of the coupled oscillators by including Gaussian 
white noise to the first variable of all the systems. In- 
terestingly, we find that the results obtained remain un- 
altered for small values of the noise intensity and the 
synchronization error, 77, follows a power law decay as 
the intensity of the noise is increased resulting in noise- 
enhanced and noise-induced synchronizations depending 
on the value of the coupling strength. It is of interest to 
note that similar noise-enhanced phase synchronization 
in two coupled noisy Rossler oscillators [llf, HI] and noise- 
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FIG. 5: (color online) Conditional Lyapunov exponents of 
the coupled (a) Rossler and, (b) Lorenz systems for different 
values of the number of oscillators in the array as a function 
of the coupling strength in the presence of Gaussian white 
noise with intensity Do = 0.001. 



induced phase synchronization in two coupled Rossler 
and Lorenz oscillators [ll|, have been observed. 

In particular, we have included the Gaussian white 
noise, y/2aDo£,{t) with a = 0.01 to the x variable of 
all the coupled systems including the drive after every 
time step, where Do is the noise intensity and £,{t) is the 
Gaussian white noise. We have calculated the CLEs for 
both the Rossler and Lorenz oscillators as in the previ- 
ous section by including a small noise with noise intensity 
Do = 0.001, which are plotted in Figs.^a) and Mjo), re- 
spectively. It is evident from this figure that the critical 
values of e for the chosen value of the noise intensity re- 
main almost the same as in Figs. [3] and |3]for the Rossler 
and Lorenz oscillators, respectively. Hence, the expo- 
nential relation between the number of oscillators in the 
ring and their corresponding Ec remains the same in the 
presence of small noise also. 

Furthermore, we have calculated the synchronization 
error rj by increasing the noise intensity Do at the thresh- 
old value of e, namely ec, shown in Table HI The average 
synchronization error (ry), where (•) denotes the time av- 
erage over 500 000 time steps, for different values of N in 
the ring as a function of noise intensity Do is shown in 
Figs, inja) and IHb) for the Rossler and Lorenz systems, 
respectively. The average synchronization error for both 
the Rossler and Lorenz systems for different numbers of 
oscillators in the ring follows a power law decay as a func- 
tion of the noise intensity Dq beyond certain threshold 
values of Do, i.e, Do > 0.01. 

In our studies, we have fixed the value of the coupling 
strength e at the critical coupling strength Sc as shown 
in Tables for different values of N. On increasing the 
noise intensity Do, the synchronization error follows a 
power-law decay as a function of Do after certain thresh- 
old value of noise intensity. This naturally corresponds to 
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FIG. 6: (color online) Time averaged synchronization error 
{■q) for different numbers of oscillators in the ring as a function 
of the noise intensity Do for a fixed value of Ec displaying a 
power law decay for (a) the Rossler oscillators in the range 
Do G (0.01, 1), and (b) the Lorenz oscillators in the range 
Do e (0.01, 10). 

a noise-enhanced synchronization and also confirms the 
robustness of the synchronous state. Interestingly, we 
have also observed that by fixing e at a value less than 
that of the synchronization threshold Sc, noise can induce 
synchronization between the ring and the drive oscillator. 
The synchronization error again shows a power-law de- 
cay as a function of the noise intensity (exactly similar 
to Figs. |6] omitted here to avoid repetition), exhibiting 
noise-induced synchronization. Thus we have observed 
both the phenomena of noise-enhanced and noise-induced 
synchronizations depending upon the choice of e. 



III. OVERCOMING SIZE INSTABILITY BY 
INTRODUCING ADDITIONAL COUPLINGS 

In the above, wc have pointed out that, when one con- 
siders an array of diffusively coupled chaotic oscillators 
(with a ring geometry) driven by an external identical 
forcing in the drive-response configuration, the critical 
coupling strength of the external drive required to syn- 
chronize the response array varies exponentially with the 
number of oscillators within the synchronization regime. 
As a consequence of the exponential relation, the number 
of oscillators in the response array that can be synchro- 
nized with the external drive is limited to 4 or 5 because 
one requires a high coupling strength, which results in 
desynchronization above a certain threshold value of the 
coupling strength due to size instability. However, we 
wish to point out here that it is possible to increase the 
number of oscillators in the response array, which are 
synchronized with the external drive by increasing the 
number of couplings. In Fig. [7] we show a schematic di- 
agram for the realization of the drive-response configu- 
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FIG. 7: (color online) Schematic diagram of the ring of diffu- 
sively coupled oscillators driven by an external identical os- 
cillator with additional coupling in (a) first neighbor and (b) 
second neighbor. 



ration vifith more than one coupling. In this section we 
explore the possibility of increasing the number of os- 
cillators that evolve in synchrony with the drive and the 
robustness of the scaling exponent 7 when one introduces 
more number of couplings between the drive and the re- 
sponse array. 

In order to study the effect of making more number of 
couplings between the drive and the response array, we 
again consider the coupled Rossler systems as discussed 
earlier in Sec. Ill Al The drive system is assumed to follow 
the same set of equations ([T]) as before. Then the gov- 
erning equations for the response array can be written as 



+ d{xj+i + Xj^i - 2xj) +Pj{xo - Xj), 



Zj = 5 + Zj (xj — c) , j 
where 



1,2, ■ 



(8a) 
(8b) 
(8c) 



is easy to see that the critical coupling strength required 
to synchronize all the three oscillators in the array gets 
reduced to « 0.311 from Sc « 0.407 (see Table 
Similarly, up to A'^ = 6 the second additional coupling 
helps to reduce the critical values of e as tabulated in 
Table |lT] for the coupled Rossler oscillators. The same 



TABLE II: Critical coupling strength Sc for second coupling 
at the first neighbour 



Rossler system 


Lorenz system 


N Ec 


iV 


£c 


3 0.311 


3 


9.189 


4 0.441 


4 


12.895 


5 0.647 


5 


18.435 


6 1.056 


6 


28.297 



analysis can be extended to the Lorenz oscillators 
d?]) as well. Again the results are tabulated in Table |lll 
It may be noted that for such a second additional cou- 
pling one additional oscillator, namely TV = 6, in the 
ring can be synchronized compared to the case of a sin- 
gle coupling (see Table |l|). Further, the critical values of 
e in Table HI] again show an exponential relation with the 
number of synchronized oscillators in the array as shown 
in Figs. ISl^a) and HJb) for the Rossler and Lorenz os- 
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Pj 



k=0 



kl + l 



1 for i = kl 
otherwise. 



1, 



(9) 



and £ corresponds to the coupling strength. Here m de- 
notes the number of couplings and I represents the l-ih 
neighbor of the first oscillator in the response array. For 
example, / = 1 for the first neighbor, I — 2 for the second 
neighbor and so on. 



A. Second coupling at the first neighbor 

First let us consider the case in which there is a sec- 
ond coupling (of the same strength as the first one) be- 
tween the external drive and the response array in the 
immediate neighborhood of the first oscillator, that is, 
Z = 1 in Eq. ([9|) , which is already coupled to the drive [cf. 
Fig.[7|a)]. For simplicity, we consider the coupled Rossler 
systems ([1]) and ([S]) with = 3 with coupling at the first 
and second oscillators in the array. Now, interestingly it 



FIG. 8: (color online) The variation of the critical coupling 
strength, Ec, as a function of the number of oscillators, N, for 
second additional coupling at the first neighbor of the first 
coupling for (a) Rossler system (eo = 0.0788 and 7 = 0.43303) 
and (b) Lorenz system (eo = 2.7351 and 7 = 0.3878). Filled 
circles corresponds to the numerical data and solid lines are 
the plot of Ec using relation ((5]). 

cillators, respectively. Here the proportionality constant 
eo = 0.0788 and the scaling exponent 7 = 0.4303 for the 
Rossler oscillators and eg = 2.7351 and 7 = 0.3878 for 
the Lorenz oscillators. 

The conditional Lyapunov exponents of the coupled 
Rossler and Lorenz oscillators for N ~ 3, 4, 5 and 6 oscil- 
lators in the ring as a function of e is shown in Figs. [9l^a) 
and[ni[b), respectively. It is to be noted that all the CLEs 
change their values from positive to negative near the 
critical values of e as in Table [TTl 

Wc have also investigated the effect of noise in the 
present case as well. The CLEs of both the coupled 
Rossler and Lorenz oscillators for different N values in 
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FIG. 9: (color online) Conditional Lyapunov exponents of the 
(a) Rossler and, (b) Lorenz systems for different values of the 
numbers of oscillators (A'') in the array as a function of the 
couphng strength with the second additional coupling at the 
first neighbor of the first couphng. 
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FIG. 10: (color online) Conditional Lyapunov exponents of 
the (a) Rossler and, (b) Lorenz systems for different values of 
the number of oscillators (A'') in the array as a function of the 
coupling strength with the second additional coupling at the 
first neighbor of the first coupling in the presence of Gaussian 
white noise with intensity Do = 0.001. 



the ring for the noise intensity Dq = 0.001 is plotted in 
Figs. [Tol a) and lTor b) as a function of e. Again the CLEs 
change their signs almost at the same critical values of 
e as the CLEs of the coupled oscillators without noise 
[Figs. [HJa) and EJb)] thereby preserving the same expo- 
nential relation between the critical values of e and the 
number of oscillators synchronized in the ring. Further, 
the synchronization error (rj) follows a power law decay as 
a function of the noise intensity Dq for both the Rossler 
and Lorenz oscillators as shown in Figs. [TIT a) andfTTTb) 
beyond certain threshold values of Dq for a fixed value 



FIG. 11: (color online) Time averaged synchronization er- 
ror (ri) for different number of oscillators in the ring as 
a function of noise intensity Dq for fixed value of Sc dis- 
playing a power-law decay (a) the Rossler oscillators in the 
range Do € (0.01, 1), and the Lorenz oscillators in the range 
Do G (0.01, 10) with the second additional coupling at the 
first neighbor of the first coupling. 

of £c given in Table |lT] illustrating noise-enhanced syn- 
chronization. Further, a similar relation can be obtained 
between (rj) and Dq for values oi e < Sc confirming the 
existence of noise-induced synchronizations. 

B. Second coupling at the second neighbor 

Next, by introducing a second coupling at the second 
neighbor (instead of the first one), ^ = 2, of the first cou- 
pling, wc find that one can synchronize two additional 
oscillators in the ring compared to the number of syn- 
chronized oscillators with a single coupling. Likewise, 
introducing the second coupling at the third neighbor 
[cf. Fig. [Zl^b)] of the first oscillator {I = 3) in the array 

TABLE III: Critical coupling strength Ec for second coupling 
at the second neighbour. 



Rossler system 


Lorenz system 


N Ec 


iV 


Ec 


4 0.405 


4 


12.212 


5 0.563 


5 


16.231 


6 0.814 


6 


21.883 


7 1.156 


7 


31.987 



results in increasing the number of oscillators that are 
synchronized with the drive by 3. Similarly, introducing 
the second coupling at the A^-th neighbor will increase 
the synchronized oscillators in the ring by N. Thus from 
the Table U for single coupling, we can realize that it is 
possible to synchronize up to 10 oscillators (instead of 
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5) with the introduction of a second couphng at reduced 
critical couphng strength. 

The critical values of the coupling strength Ec and their 
corresponding number of synchronized oscillators in the 
ring for both the Rossler and Lorenz oscillators for the 
second coupling at the second neighbor are tabulated in 
Table IIIIl which again displays an exponential relation 
as shown in Figs. [T^a) and [T^b). The proportional- 
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FIG. 12: (color online) The variation of the critical coupling 
strength, Ec, as a function of the number of oscillators, A'^, for 
second additional coupling at the second neighbor of the first 
coupling for (a) Rossler system (eo = 0.0972 and 7 = 0.3536) 
and (b) Lorenz system (eo = 3.0908 and 7 = 0.3321). Filled 
circles corresponds to the numerical data and solid lines are 
the plot of £c using relation ([5]). 

ity constant is estimated as Eq = 0.0972 and the scal- 
ing exponent 7 = 0.3536 for the Rossler oscillators and 
£0 = 3.0908 and 7 — 0.3321 for the Lorenz oscillators. 
These critical values of e are also confirmed by the tran- 
sitions in the value of the CLEs of the coupled Rossler 
and Lorenz oscillators as shown in Figs. ITST a) and lTST b). 
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FIG. 13: (color online) Conditional Lyapunov exponents of 
the (a) Rossler and, (b) Lorenz systems for different values 
of the number of oscillators in the array as a function of the 
coupling strength with the second additional coupling at the 
second neighbor of the first coupling. 
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FIG. 14: (color online) Conditional Lyapunov exponents of 
the (a) Rossler and, (b) Lorenz systems for different values 
of the number of oscillators in the array as a function of the 
coupling strength with the second additional coupling at the 
second neighbor of the first coupling in the presence of the 
Gaussian white noise with noise intensity Do = 0.001. 



0.001 for N ~ 4, 5, 6 and 7 oscillators in the ring are 
shown in Figs. E^a) and [TW b) . It is evident from this 
figure that the critical value of e for small noise intensity 
is almost the same as that of the oscillators without noise 
[cf. Figs. [T^Ta) and [T5Fb')]. indicating the robustness of 
the synchronous state and the exponential relation ([5]) 
with the addition of small noise. Further the power law 




FIG. 15: (color online) Time averaged synchronization error 
(7;) for different number of oscillators in the ring as a function 
of noise intensity Do for fixed value of Ec displaying power- 
law decay for (a) the Rossler oscillators in the range Do £ 
(0.01, 1), and (b) the Lorenz oscillators in the range Do £ 
(0.01, 10) with the second additional coupling at the second 
neighbor of the first coupling. 



The conditional Lyapunov exponents of both the os- 
cillators as a function of e with the noise intensity Dq = 



decay of the synchronization error as a function of the 
noise intensity [Figs. [T5lja)-(b)] beyond certain threshold 
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values of Dq for the value of Ec for both the oscillators 
confirming the existence of noise-enhanced synchroniza- 
tion. Exactly similar figures can be obtained by fixing 
the coupling strength to be less than Ec indicating the 
existence of noise-induced synchronization. 



C. Couplings at the A^-th neighbors 

Let Nc be the maximum number of oscillators in the 
response array that are synchronized with the external 
drive for a given Ec for m ~ I. Then, by introducing 
a second coupling, (m = 2) in eq. (|9]), at Nc + I one 
can synchronize up to a maximum of 2Nc oscillators in 
the response array with the external drive for the same 
coupling strength Ec- One can also introduce a third 
coupling in the response array, (to = 3) in eq. in the 
same fashion in which the second coupling is introduced 
where one can synchronize a maximum 3Nc oscillators 
in the response array. In this way one can increase the 
number of oscillators (size) in the response array that 
evolve in synchrony with the external drive to mNc by 
introducing to couplings at the oscillator index 1, Nc + l, 
2Nc + 1, . . . , (to — l)A'c + 1, respectively, with the same 
critical coupling e^- 

In the following wc shall illustrate this result using the 
second and third couplings at the oscillators with the 
indices Nc + l and 2Nc + 1 , respectively. It is known from 
Section |TT] that the maximum number of oscillators that 
can be synchronized with the drive with a single coupling 
is A^c = 5. The CLEs of the Rossler and Lorcnz oscillators 
with the second coupling, (to = 2), at the oscillator in 
the ring with the index iVc-l-l = 6 and the third coupling, 
(to = 3), at the oscillator index 2Nc + 1 = 11 along with 
the CLE of the Nc oscillators in the ring with a single 
coupling are shown in Figs. lTCT a) and ll6f b) as a function 
of e. It clearly shows that the mNc oscillators in the 
ring with couplings at to = 1,2 and 3 are synchronized 
exactly at the same critical value of the coupling strength 
Ec- Hence, it is evident that the number of synchronized 
oscillators can be increased beyond the size instability to 
any desired amount by introducing additional couplings 
at the {mNc + l)-th oscillators for the same value of Ec- 

We have also examined the effect of noise as in the 
previous cases. The CLEs of the Rossler and Lorenz 
oscillators with the second and the third couplings at 
the oscillators with the indices Nc + l and 2Nc + 1 with 
the noise intensity Dq = 0.001 for fixed value of Ec are 
shown in Figs. [TWa) and [TBTb) . respectively. The syn- 
chronization error also displays a power law decay as a 
function of the noise intensity Dq as shown in Figs.fTWa') 
and [TW b) beyond after certain threshold values of Dq 
for both the coupled Rossler and Lorenz oscillators con- 
firming the robustness of the synchronous states and the 
noise-enhanced synchronization. Further these systems 
of coupled oscillators also display noise-induced synchro- 
nization by exhibiting similar figures as a function of Dq 
for E less than Ec- 
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FIG. 16: (color online) Conditional Lyapunov exponents of 
the (a) Rossler and, (b) Lorenz system for different values 
of the number of oscillators in the array as a function of the 
coupling strength with the second and the third additional 
couplings at the A'^-th and 2A''-th neighbor, respectively, of 
the first coupling. 
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FIG. 17: (color online) Conditional Lyapunov exponents of 
the (a) Rossler and, (b) Lorenz system for different values 
of the number of oscillators in the array as a function of the 
coupling strength with second additional coupling at the Ai'-th 
and 2A''-th neighbors, respectively, of the first coupling in the 
presence of Gaussian white noise with intensity Do — 0.001. 



D. Effect of additional couplings 

In addition to the effect of increasing the number of 
synchronized oscillators in the array beyond the size 
instability limit by introducing additional couplings at 
{mNc + l)-th oscillators for to = 1, 2, • • • it reduces the 
required coupling strength exponentially to synchronize 
the fixed number of oscillators in the array. For instance, 
we have fixed the number of oscillators in the array to 
be A'^ = Nc = 5 and increase the number of couplings 
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between the drive and the array from iV = 1 to iV = 5-th 
oscillator and estimate the critical coupling strength Sc 
required to synchronize the N ~ Nc ^ 5 oscillators in 
the array for each additional coupling. The estimated £c 
is plotted as a function of the number of additional cou- 
plings in the array in Figs. [18] for both the Rossler and 
the Lorenz systems, which establishes an exponential de- 
crease of the required coupling strength to synchronize 
the fixed number of oscillators in the array by introduc- 
ing additional couplings between the drive and the array. 
It may be noted that the critical coupling strength follows 
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FIG. 18: (color online) Critical coupling strength as a function 
of the number of additional couplings between the drive and 
the response array for a fixed number of oscillators in the 
array displaying an exponential decay of Sc as a function of 
number of additional couplings (a) Rossler and, (b) Lorenz 
system. 

an exponential relation with the number of couplings as 
Ec ~ eo6xp(— 7 m), with m being the number of cou- 
plings. The constants, in the case of Rossler equations, 
turn out to be eo = 4.4746 and 7 = 0.8819 while they 
are £0 = 105.6375 and 7 = 0.7069 for Lorenz equations. 



IV. SUMMARY AND CONCLUSION 

In this paper, we have studied chaos synchronization 
in arrays of diffusively coupled nonlinear oscillators with 
a ring geometry driven externally by an identical oscilla- 
tor. In particular, we have shown that the critical cou- 
pling strength required to synchronize the array with the 
external drive increases exponentially with a scaling ex- 
ponent 7 G (0.3, 0.5) as a function of the number Nc of 
the oscillators in the array. We have pointed out that 
as a consequence of the exponential relation, the maxi- 
mum number of oscillators in the array that can evolve 
in synchrony with the external drive is limited. Further, 
we have shown that by introducing additional couplings 
between the external drive and the array at {mNc + l)-th 
oscillators in the ring, one can proportionately increase 
the maximum number of oscillators that can evolve in 
synchrony with the drive. Further, we have obtained the 
same exponential relation connecting the critical coupling 
strength and the number of oscillators even after intro- 
ducing the additional number of couplings. Furthermore, 
we have found that Sc establishes an exponential decay as 
a function of the number of additional couplings between 
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FIG. 19: (color online) Time averaged synchronization error 
(7;) for different numbers of oscillators in the ring as a function 
of noise intensity Dg for the fixed value of Sc displaying power- 
law decay for (a) the Rossler oscillators in the range Do £ 
(0.01, 1), and (b) the Lorenz oscillators in the range Do £ 
(0.01, 10) with second additional coupling at the A^-th and 
2A^-th neighbor, respectively, of the first coupling. 



the drive and the response array for a fixed number of 
oscillators in the array. We have also examined the ro- 
bustness of the results against noise of small intensity and 
found that the synchronization error displays a power law 
decay as a function of the noise intensity at e = £c indi- 
cating the existence of noise-enhanced and noise-induced 
synchronization for £ < £c in all the cases. 

In addition, we have also obtained similar results as 
above in other ubiquitous coupled nonlinear oscillators 
such as coupled MLC circuits, Chua's circuits and Sprott 
oscillators as well, and also in a discrete system, namely 
coupled logistic maps, thus confirming the universality of 
the above results. One can also extend the same type of 
analysis with other coupling configurations such as star- 
type, unidirectional, global, weighted coupling configura- 
tions, 2-d, 3-d lattices, etc. We believe that our results 
shed more light on controllability and synchronizability 
of networks by introducing additional couplings at ap- 
propriate oscillators/nodes with the less cost in terms 
of the coupling strength. Further investigations can be 
extended to networks, in particular to network with com- 
munity structure using pinning control and also with de- 
lay coupling. 
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